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Abstract 

B meson rare decays (B — ► K(K*)ll and B — > K*j) are analyzed in the 
framework of effective field theory of heavy quarks. The semileptonic and 
penguin type form factors for these decays are calculated by using the light 
cone sum rules method at the leading order of 1/rriQ expansion. Four exact 
relations between the two types of form factors are obtained at the leading 
order of 1/rriQ expansion. Of particular, the relations are found to hold for 
whole momentum transfer region. We also investigate the validity of the 
relations resulted from the large energy effective theory based on the general 
relations obtained in the present approach. The branching ratios of the rare 
decays are presented and their potential importance for extracting the CKM 
matrix elements and probing new physics is emphasized. 

PACS numbers: 11.20.Hv, ll.55.Hx, 12.39.Hg, 13.20.Fc, 13.20.He 
Keywords: rare decay, CKM matrix elements, form factors, effective field 
theory, light cone sum rules 
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I. INTRODUCTION 



B meson rare decays B — > K*j and B — ► K(K*)U are induced by transitions 6 — > 57 
and 6 — > sZZ via penguin loop diagrams in the quark level, which are usually called flavor- 
changing-neutral-current (FCNC) processes. In the Standard Model (SM), B meson rare 
decays may provide a quantitative way to determine the CKM matrix elements Vu, Vt s and 
V t b- In the SM, FCNC transitions are forbidden at tree level. They can be induced only 
starting at 1-loop order, which makes their rates for decaying be sensitive to probe new 
physics. For these reasons, they have long been hot subjects in both experimental and 
theoretical studies. 

Nevertheless, one remains facing difficulties in studying exclusive B meson rare decays 
due to the requirment of explicit calculations for the relevant form factors which involve in 
long distance ingredients that can not be calculated via QCD perturbative theory. Some 
reasonable nonperturbative methods, such as QCD sum rules, lattice simulations and phe- 
nomenological models, have been developped to estimate the long distance effects. 

The rare decays have been studied by using light cone sum rules (LCSR) in full QCD 
theory P],0,]3|,0J] . As B meson can be treated as a heavy meson containing a single heavy 
quark, it is of interest to apply for the heavy quark effective field theory (HQEFT) to deal 
with the B meson rare decays. In this paper, we shall use the framework and normalization 
derived in Refs. Hm^] , which has been applied in Refs. fl|TCi|| to investigate the exclusive 
semileptionic B decays into ir and p by using the LCSR method and obtained quite rea- 
sonable results. A more general study on exclusive semileptonic decays of heavy to light 
mesons within the framework of HQEFT has recently been carried out in Note that 
as we only keep to the leading order contributions in the expansion of 1/mg, the results 
and conclusions are actually independent of any framework of effective field theory of heavy 
quarks. 

We shall focus on in this paper the LCSR calculations of the form factors for the exclusive 
B meson rare decays B — > K*^ and B — > K(K*)ll within the framework of HQEFT. In 
section |H], we first present the hadronic matrix elements in the framework of HQEFT and 
derive a set of formulae for the relevant form factors. In section [TTT| , LCSR approach is 
applied to the relevant correlator functions in HQEFT. We then obtain four interesting 
relations among semileptonic type form factors and penguin type ones. Of interest, these 
relations are found to hold in whole momentum transfer region in the infinite mass limit of 
heavy quark or at the leading order of 1/w,q expansion in HQEFT. Obviously, at zero recoil 
of the final light meson, these relations recover the so-called Isgur-Wise relations [fT^] . In 
early time, the Isgur-Wise relations were conjectured to be also valid in the region of large 
recoil in ref. [O. Late on, these relations were really shown in the quark model to 



hold at large recoil. In particular, one of the relations concerning the form factor in radiative 
decays was shown to hold in the whole momentum transfer by using QCD sum rule approach 
[16| . Numerical analysis of form factors is presented in section |IV|. Recently, one developed 



the so-called large energy effective theory (LEET) in which more relations were obtained 
near large recoil due to additional symmetries in large energy limit while some of 

the relations were found to be broken down by QCD corrections ||19|| . Since all relations in 
Ref. |I7j were put forward following from the combination of heavy quark effective theory 
(HQET) and LEET and moreover LEET is compatible with LCSR, our present method 
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provides a reliable and important way to check the validity of LEET relations. So a detailed 
discussion and comparison will be presented on the basis of our HQEFT calaulation in 
section [V]. In section |VT], we give the relevant branching ratios for the B meson rare decays 
B — > K*"f and B — > K(K*)ll. A brief summary is outlined in section [VII. 



II. GENERAL DESCRIPTION OF MATRIX ELEMENTS IN HQEFT 

The transition matrix elements responsible for the B meson rare decays B — > K*j 
and B — > K(K*)ll may be grouped into two types: semileptonic and penguin ones. The 
semileptonic ones are defined as 

(K(p)\s^b\B(p + q)) = 2/ + (gV + (U(q 2 ) + f-(q 2 )W (2-1) 

for B to K decays and 

(K*(p, e*)\sY(l - j 5 )b\B(p + q)) = -%{m B + m^AM 2 )^ 

^ ^ (e* ■ (p + q)){2p + qT + i A3 }_ q2) {e* ■ (p + q))<f 
m B + m K * m B + m K * 

+ e^elip + q) pPl (2.2) 

for B to K* decays. In this paper, we take e i23 = 1 and 7 5 = 75 = Z7°7 1 7 2 7 3 . 
For convenience, we may define a form factor Ao(q 2 ) as 



a ( 2n 2{m B + m K *)m K * - 2 2 
t / 2\ ( m B + «^*)A(<? 2 ) - (m B - m K *)v4 2 (g 2 ) 

A 3W J — ~ 

2m K * 

A Q (0) = A 3 (0) (2.3) 

A (q 2 ) will directly enter into contributions to the relevant branching ratios. 
The penguin matrix elements can be written as 

(K(p)\aa^q v (l + l 5 )b\B(p + q)) = 1 ^ {q 2 (2p + qf - (m| - m 2 K )q»} (2.4) 

for B to K decays and 

(^(p,e*)|sa^(l + 7 5 )&|5(p + g)) = -ie^e;(p + g)^ 7 2T 1 (g 2 ) 
+T 2 (g 2 ){(m| - m 2 K *y» - (e* ■ (p + q))(2p + gf} 

+T 3 (q 2 )(e* ■ (p + q)){q» - 2 ^ 2 (2p + (2.5) 

— * 

for B to if* decays. 

In the above definitions, p is the momentum of the light meson K or K* . e* is the 
polarization vector of K* meson, and q is the momentum transfer. f± and fx are the B 
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to K semileptonic and penguin transition form factors respectively. Aj(i=0,l,2), V and 
Tj(i=l,2,3) are the corresponding ones for B to K* transitions. 

To be convenient for making Borel transformation which helps to suppress the contribu- 
tions from the possible higher states of bottom mesons, we may change (2.4) and (2.5) into 
the following forms 

(K(p)\sa^p v (l + ^)b\B{p + q)) = i /T ^ 2) {(q ■ p)(2p + q)» - ((2p + q) ■ p)(f } (2.6) 



(K*(p,e*)\sa^p„(l + ^)b\B(p + q)) 
—ie 



"'*V (p + Q)0Py{ m|+ ^ + g2 T 1 (g 2 ) + ™k. 



m 



K 



: T 2 (g 2 )} 



+e^{(q-p) 



rn 



B 



m 



K 



-T 2 (q 2 )-[(q-p) 



m B — m 



K* 



{2p + q)-p)T 1 (q 2 )} 



+q h 



(e*-(p + q))((2p + q)-p) 



m 



B 



m 



W) + 



m, 



niK *-{T 2 {q*)-T x {q*))} 

-W)-T^))} (2.7) 



K* 



_ (2p + 9 ),(lM|±fc) { T 3 tf) 



m B — m 



m B — m 



Is 



When applying for the HQEFT to evaluate the matrix elements, they can be expanded 
into the powers of 1/ mq and also be simply expressed by a set of heavy spin- flavor indepen- 
dent universal wave functions ||j|||,|10|] . It is convenient to adopt the following normalization 
which relates matrix elements in full QCD with the ones in HQEFT || 



[K\sTb\B) 



-.{(K\sTb v \B v ) + 0(l/m b )} 



(2.8) 



where k represents K(p) or K*(p,e*). The notation b v is the effective bottom quark field. 
And A B = itlb — vn\j is the binding energy. From heavy quark symmetry, one can obtain the 



following relations |, 10, 20, 21, 22 



with 



(K(p)\sTb v \B v ) = -Tr[k{v,p)TM v ] 
(K*(p,e*)\sTb v \B v ) = -iTr[Q(v,p)TM v 



k(v,p) = rf[A(v ■ p, fi) + f>B(v -p, fi)] 

n(v,p) = Lt{v • p)f + L 2 (v -p)(v e*) + [L 3 {v ■ p)f + L A {v ■ p){v ■ e*)]f> 



(2.9) 
(2.10) 



(2.11) 
(2.12) 



and 



pf 



v ■ p 



M v = -VI 



(2.13) 
(2-14) 
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Where A, B and Li(i = 1,2,3,4) are the leading order wave functions characterizing the 
heavy to light transition matrix elements in the effective field theory. M. v is the spin wave 
function associated with the heavy meson state. The vector v M is the four-velocity of B 
meson satisfying v 2 — 1, and A is the heavy flavor independent binding energy 

A = lim A B 

mQ— +00 

which reflects only the effects arising from the light degrees of freedom in the heavy B meson. 
With eqs.(2.1-2.14), one arrives at the following expressions for the form factors 



fr(q' 

Mq' 

V(q- 

T2(q' 

T 3 (q' 
with 



1 

m B 



m B + m K 

m B ^ Ai 

2 m B A 
m B + m k * ^ A B 

2(m B + m K *) 



m B A B'(v -p) 



v ■ p 



+ 



{L^v p) + L 3 (v p)} + 



m B A ( L 2 (v -p) L 3 (v-p) - L 4 (v-p) _ 



2m% 



2m B (v ■ p) 



2(m B + m K *)> 



m B A ( L 2 (v -p) _ L 3 (v-p) - L 4 (v ■ p) _ 
1 2ml °~ '" ^ 1 



2m B (v ■ p) 



m B A m B + m K * 
\ A B m B {v-p) 



L 3 (vp) + 



m B A L[(vp) L' 3 (vp) 

\ I 1 J I 

\ A B m B vp 



= 2, 



m B A 



\ A B m\- m\. 



■{{m B - v •p)L' 1 (vp) + 



m B v ■ p - m K » 
v ■ p 



L' 3 (vp)} + 



m B A^ L[(v-p) + L' 3 (v-p) 



m 



m 



m B 



v ■ p 



y = vp 



m B v ■ p 



^L' 4 (vp)} + 



m\ + m 2 K — q 2 
2m r 



(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 
(2.20) 
(2.21) 
(2.22) 
(2.23) 



(2.24) 



denoting the energy of the final light meson. 

In the above formulae the dots denote possible higher order I/tuq contributions that are 
neglected in this paper. B'{v ■ p) and L'^v ■ p) are in general different from B(v • p) and 
Li(v • p)(i = 1, 2, 3, 4) as they arise from different matrix elements. 
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III. LIGHT-CONE SUM RULES IN HQEFT 



In order to calculate the relevant hadronic matrix elements which contain nonperturba- 
tive contributions and thus make QCD perturbative method lose its power, we shall apply 
for LCSR approach. In LCSR calculation, the relevant correlation functions are expanded 
near the light cone. The light cone distribution functions are introduced to describe the 
nonperturbative effects. In searching for reasonable and stable results, the quark-hadron 
duality and Borel transformation are generally adopted (for a detailed review, one may find 
in Refs. 

The theoretical calculations can often show a simpler process in the framework of HQEFT 
than in QCD, which can explicitly be seen in Refs. [§,|TIJ where the semileptonic form factors 
for B — ► 7i and p have been evaluated. We may directly adopt the formulae in p,|TU,n| to 
B — > K transitions by simply changing the relevant quantities corresponding to the K meson 



f. 



K 



AFy Jo 



*M - "if 



dse T 



dse t 



ld_ 

y du 



92(11) - /i K (f) p (u) 



6 du 



u=l- 



■2y 



-(pK{U) + -o a-offlW j7-j2 « + ~~a ~a rff(^) 

y z ou z y A au by ou 



(3.1) 
(3.2) 



u=l-f- 



As for the B to vector Kaon meson decays, we may just adopt the definitions of K* 
meson distribution functions given in [2^ . 



(iT(p,e*M-xK^(x)|0> = -ifU(t> 



[iris ^ v PfJ,) 

N e* • x 



due-^UAu) + 



m K *x 



A T {u)) 



(p ■ x) 



-m 



K* 



1 due-^ x B T {u) + l -{tlx v 



m 



K" 



p • X JO 



(K*(p,e*)\s(-x) lfl d(x)\0) = fK*m K *iP, 



due-^ p - x C T {u)} 



x 



+ e 



A' 



p ■ x Jo 
e ■ x 



p ■ X Jo 

1 e* ■ x 
-x,,- —m 



2 ^{p-xf 



K* 



due- iip - x {4>\\{u) 
due- l ^ x g { l\u) 
due-^ p - x C(u)} 



2 2 

m K *x 



A{u)) 



o 



(K*(p,e*)\s(-x)w 5 d(x)\0) = hf K * - f^ ms + md )m K *e^ a pe*»p a xP 

-1 



(3.3) 



(3.4) 

due^ p - x g ( f{u) (3.5) 



(K*(p,e*)\s(-x)d(x)\0)=i(f^-f K 



m s + m d , 
m K * 



)m 



K* 



x) 



due 



h\>(u) 



(3.6) 



with £ = 2u — 1. In the above definitions 4>±(u) and 4>\\{u) give the leading twist two 
distributions in the fraction of total momentum carried by quarks in the transversely and 
longitudinally polarized meson respectively. The functions g±(u) and g^\u), which de- 
scribe transverse polarizations of quarks in the longitudinally polarized meson, provide con- 
tributions of twist three, hi?'' is also twist three function. A(u), A T (u) and C(u) are twist 
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four, while Bt{u) and Ct{u) are mixing twist functions. In this paper, we only consider 
distribution functions from two-particle contributions. 

Applying the same calculation steps as in [10], we arrive at the following results 



Li(y) 



l 

4F 
1 
4 



e 2 ~^ T rdse-" T - \f K .m K *g { l\u) 

y L 



® (<*)/ \ i fx* 7712 } 
mj{* ou 2y 



+- a Uk* - Jr* ) m K*—g±[u) + — C T (u) 



(3.7) 



2y 



L 2 {y) 



e 2A */ T / dse- s ' T 



1 



fjc.m 2 K . ( l 

y 2 l 2 



, ( m K . m s + m d d 

2y 3 2y 2 m K * du 



h { {\u)) 



(3.8) 



23/ 



1 



L ^ = JF e2AB/T i dS6 ~ S/T 
, s m 2 K , d 2 



1 , , fT m s + m d d (o) 

4y 7T7.K* C'M 



^4(2/) 



1 



M - 



16y 2 9m 2 



A T («)) 



(3.9) 



u=i- 



2.U 



4F 

f K *m K » 



-1 
2/ 



A'* 



1 



1 d 



m 2 K » d 2 



K u)-g ± '{u)--—gy{u) 







IQy 2 du 2 



A{u)) 



BT(u) + -jl*(m s + m d )j^-g { l\u) 



(3.10) 



14=1- 



-'!/ 



For the penguin matrix elements, we begin with considering the following correlations 
P»( Py q) = i J d 4 xe^(K(p)|T{s(x)^V(l + 7 5 )6(x),6(0)z 7 5 rf(0)}|0) (3.11) 
V£(p,q) = i [ d 4 xe- ipB - x (K*(p,e*)\T{s(0)a^p v (l + 7 5 )6(0), 6(x)i 7 5 rf(x)}|0) (3.12) 



where Pb = P + Q- 

Making the procedures similar to Refs. PJlQl, one can easily derive the explicit expres- 
sions for B'(y) and L\{y\ which are exactly the same as B(y) and Li(y) in (3.2) and (3.7-3.10) 
at the leading order of 1/mg expansion. From (2.15) to (2.23), we obtain the following four 
interesting relations among the semileptonic type and penguin type form factors. 



f T (q 2 ) 

T2(q 2 ) 
T 3 (q 2 ) 



U+{q 2 ) - Uq 2 )) 



IrriB 

m 2 B — m 2 K * + q 2 V(q 2 ) 



+ m B +m K * Mq2] 



2ms tub + rriK* Ivtib 
2 (m B -y)(m B + m K 



m B {y 2 - m 2 K „). 



m B — m 
m 2 B — m 
m B q 2 



K* 



m B + m K * 



(3.13) 
(3.14) 
(3.15) 



A'* 



m K *A (q 2 



m 2 B — m 2 K * + q 2 
2rriBq 2 



[{m B + mK^A^q 2 ) - (m B - m K *)A 2 (q 2 )} 
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,n| + 3,<.-^ 
2m B {m B + m/<») 

Note that these relations hold in the whole momentum transfer region. In fact, they can 
simply be obtained from eqs.(2.9) and (2.10) by assuming that the light meson wave func- 
tions k(v,p) and Q(v,p) are universal for different types of the currents, namely they are 
independent of the choice of T in the currents. The second relation was actually observed 
in ref. |16|] by using QCD sum rule approach. 



In the heavy quark limit near zero recoil point (q 2 — > q^ax), the above relations were 
known as Isgur-Wise relations [O. In Refs. |2"7|,|2"8"[| , some detailed discussions were made to 



explore the Isgur-Wise relations in both small and large recoil conditions by using full LCSR 
method, where it was shown that the Isgur-Wise relations are satisfied very well at q 2 — > 
(large recoil) and hold with about 80% accuracy at large q 2 . In Ref. p9fl , conclusions were 



made that the Isgur-Wise relations are valid up to 70% in the whole q 2 region by applying 
for the three point QCD sum rules method. In the quark model, the authors of Refs. flT^Jl"5| 
concluded that Isgur-Wise relations also hold at large recoil. 

In here we have further confirmed the existed analyzes in the literature and arrived at a 
general conclusion that the relations (3.13-3.16) hold for whole allowed region of momentum 
transfer q 2 at the leading order of 1/mg expansion, which has nothing to do with the 
energy of the final light meson. This may provide a more general proof in support of the 



hypothesis made in the early time by the authors in Ref. |L3| . As a consequence, it becomes 
remarkable that one can directly read off the penguin form factors fa and Tj(i = 1,2,3) 
from the relations in eqs. (3.13-3.16) without the tedious calculations from LCSR to Borel 



transformation. In Ref. [pTTJ] , we have shown that when applying the HQEFT to heavy to 
light semileptonic decays, heavy quark expansion and heavy quark symmetry enable us to 
relate various decay channels, consequently, the theoretical analysis is much simplified and 
the number of independent functions is greatly abated though the number of independent 
functions in a single decay channel does not decrease. Here with the four exact relations, 
the number of independent variables among the form factors f± and J't is straightforwardly 
reduced to two. The number of independent variables in the rare decay B — > K*ll is then 
reduced from seven form factors to four form factors. One sees that for the B meson rare 
decays the heavy quark expansion exhibits its more powerful advantages, the number of 
independent functions is found to be largely reduced even in a single decay channel. 



IV. NUMERICAL ANALYSIS AND RESULTS FOR FORM FACTORS 



The light cone wave functions play an important role for a precise calculation of form fac- 
tors. They have been studied by several groups. We shall use the results given in P^J5D|j5T]J5^ 



for K and the ones in p7| for K* meson in our following analyses. The asymptotic form 
and the scale dependence of them are taken from perturbative QCD calculations given in 

We shall first present all two-particle light cone amplitudes of Kaon appearing in (3.1) 
and (3.2) 



>K 



[u, 11) = 6u(l — u) 



a 1 (^)(3(2«-l)) 



15 

a2(^)(y(2w 



3, 
2' 



S 



35 1 ^ I*) 

+o 8 (A*)(y (2u - l) 3 - y (2u - 1)) + a 4 (/i)^(21(2 M - l) 4 - U(2u - l) 2 + 1) 

b p (u, fi) = l + ~B 2 (/i) [3(2m - l) 2 - 1] + ~B 4 (/i) [35(2m - l) 4 - 30(2w - l) 2 + 3] 

1 8 

3 1 F5 

^(u,/x) = 6u(l - + -<7 2 (/i)[5(2u - l) 2 - 1] + — C 4 (/x)[21(2u - l) 4 - 14(2u - l) 2 + 1]} 

Z o 



^(u, M ) = -S 2 (^)u 2 (l - uf + ^e(/x)5 2 (/x) 



6 

u(l - u)(2 + 13m(1 -u) + 10M 3 (logw)(2 - 3u + -u 

5 



+ 10(1 - uf i (log(l - u))(2 - 3(1 - «) + |(1 - u) 2 )) 

5 

^ 2 K / i) = H ( 5 2 ( /U ) M (l- M )(2 M -l) (4 
For the relevant parameters, we will take the following values in our numerical calculations 

ai{p>b) = 0.15, a 2 (fib) = 0.16, a 3 (/i b ) = 0.05, a 4 (/x&) = 0.06, 
5 2 (/i fe ) = 0.29, 5 4 (/i 6 ) = 0.58, C 2 (a* 6 ) = 0.059, C 4 (/i b ) = 0.034, 

5 2 (/i 6 ) = 0.17GeV 2 , e(fi b ) = 0.36 (4.2) 



where we choose fXb = \Jrn 2 B ~ m l ~ 2.4GeV which is the appropriate scale characterizing 
the typical virtuality of the b quark. The SU(3) flavor violation effects of Kaon has been 
taken into account in 4>k(u, /x) with non- vanishing coefficients ai(/ife) and as(fib)- The SU(3) 
flavor breaking effects in the higher twist amplitudes are neglected in this paper as they only 



have a small contribution ||3"5| 



For the light vector K* meson wave functions, we use the expressions given in |27 
The light meson SU(3) flavor breaking effects are taken into account in the leading twist 
distributions and partially in the twist three, but ignored in twist four. 

Other parameters needed are listed in the following : 

m B = 5.28GeV, m h = 4.75GeV, m K = 0.49GeV, m K * = 0.89GeV, 
A B = 0.53GeV, A = (0.53 ± 0.08)GeV, F = (0.30 ± 0.06)GeV 3/2 , 
f K = 0.16GeV, f K . = (226 ± 28)MeV, fa(ji b ) = (175 ± 9)MeV, 

HkM = « 2 ' 02GeV ( 4 - 3 ) 

As for the parameters So and T, according to LCSR criterion that the contributions from 
both the higher states and higher twist four distribution functions should not be larger than 
30%. Thus we choose the region for T to be lGeV < T < 3GeV. In this region the curves 
of f± and fx for B to K transitions and A^i = 0, 1, 2), V and T;(z = 1, 2, 3) for B to K* 
transitions become most stable at the threshold energy so = 2.1 — 2.7(2.4 ± 0.3)GeV and 
So = 1-8 — 2.4(2.1 ±0.3)GeV respectively, which may be seen explicitly from Fig.l to Fig. 10. 

It is known that the light cone expansion and the sum rule method may be broken down 



at large momentum transfer (practically as q approaches near half of m 2 ) p4"| , which may 
be seen explicitly from Fig. 11-20. Where the curves of form factors calculated from LCSR 
likely become unstable at large q 2 region. Thus for the behavior of the form factors in 
the whole kinematically accessible region, we may use the following parametrization for the 
transfer momentum dependence of the form factors 
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where F(q 2 ) can be any of the form factors /+, /_, fx, A^i = 0, 1, 2), V and T^i = 1, 2, 3). 
We directly use LCSR predictions to fit the parameters. This is because K* and K mesons 
contain a relative heavy strange quark (in comparison with the u and d quarks), so that 
they are comparatively heavy and shorten the kinematically allowed ranges of B to K* and 
K decays in comparison with the ranges for B to 7r decays. Therefore the sum rules are 
expected to yield reasonable values for most allowed region of q 2 in the K* and K meson 
cases. 

With the above analyses, we are then able to fix the three parameters for each form 
factor. We plot in Fig. 11-20 the form factors as functions of q 2 with different threshold 
energy s = 2.1, 2.4, 2.7GeV for B to K transitions and s = 1.8, 2.1, 2.4GeV for B to 
K* transitions at T = 2.0GeV. Our numerical results of form factors at q 2 = are given 
in Table 1, where the uncertainties of the form factors arise from the uncertainties of the 
threshold energy so, the Borel parameter T and the parameters in (4.3). In our values, the 
uncertainties coming from so are about 10%, from T and the parameters in (4.3) are around 
15%. The total uncertainties are up to 25%. For comparison, we list in table 1 the numerical 
results obtained from other approaches: QCD LCSR, LEET, quark model (QM), lattice, 
three point sum rules (SR) and PQCD calculation. 

Table 1. Values for form factors of B -> Kll,B -> K*ll and B -> K*j 
at q 2 = 0. 





present 


LCSR |I] 


LEET [0 


QM 


lattice PB| 


lattice PI 


SR |2| 


PQCD 


/+(«) 


U.404 _ .075 


0.319 tZt 




0.36 




0.30(4) 


0.25 




/o(0) 


U.^O^ _ .075 


0.319 t^i 




0.36 




0.30(4) 


0.25 




Mo) 


0-447 2S8S 


n o« +0.016 
U.OdO -0.055 




0.35 




0.29(6) 


0.14 




A (0) 


0.468 


n 471 +u.i'27 

U-^' 1 -0.059 




0.45 


0.32 




0.30 


0.407 


Mo) 


n op-n +U.U68 
U.-30U _o.089 


n QQ7 +0.048 
U.-30^ _o.043 


0.27 ±0.03 


0.36 


0.28 




0.37 


0.266 


A 2 (0) 


U.OUZ _ .080 


n 9 89 +0.038 
U.ZOZ -Q.036 




0.32 






0.40 


0.202 


V(0) 


426 +u - uvu 

U.^ZO _ .098 


0.457 ™ 


0.36 ±0.04 


0.44 


0.38 




0.47 


0.355 


li(0) 


U.OOZ -0.093 


0.379 ^ 


0.31 ±0.02 


0.39 


0.32 




0.38 


0.315 


r 2 (o) 


n ^89 +U.068 
U.OOZ _ .093 


n ^7Q + U - U58 

U.O / J _o.045 




0.39 


0.32 




0.38 


0.315 


r s (o) 


o 9m +U - U45 

U.ZOO _ .063 


U.ZOU _o.016 




0.27 






1.4 


0.207 



V. ON VALIDITY OF RELATIONS IN LEET 



Recently, it was noted in Ref. |T7] that by applying HQET to the initial heavy meson and 
meanwhile adopting LEET to describe the final light meson, more relations were obtained 
at the leading order of heavy quark mass and large energy expansion. Eventually, the total 
number of independent form factors is reduced to three near the large recoil point. While 



those symmetry relations were shown |y| to be broken down when radiative corrections are 
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considered. The contributions of the second order in the ratio of the light meson mass to 
the large recoil energy were taken into account in Ref. fll8f| . 

To compare our results with those in Ref. [17] and also to see how much accuracy of the 
large energy effective theory would be, we first introduce the following definitions 



C(m B ,y 
(A(m B ,y 
C\\(mB,y 
(±(m B ,y 
(i(m B ,y 
(2(m B ,y 



f±(q 2 ) 
fo(q 2 ) 

f T (q 2 ) 
A (q 2 ) 
Ai(q 2 ) 



A,(a 2 ) 



V(q 2 ) 
Ti(q 2 ) 

T2(q 2 ) 



(A(m B , y)±( 



m^A 1 



b y 



B{y) 



mgA 1 



b m B 



My) 



m B A 1 



B m K * 



My) - 



y 



m B m K * 



My)) 



m B A L 3 (y) 



b y 



m B ALi(y) 



b y 



m B A y 



My) 



\J A B m B m K » 

With these functions, the form factors in eqs. (2. 15-2.23) can be expressed as 

(m B ,y) 



f-{q 2 ) + U(q 2 ) 



'I - 



- — 2~)(( m B, y) + (l + —2-^ 



m B — m 



K 



m B — m 



-)(A(m B ,y) 



K 



1 + —)((m B} y) 



(1 



m B 

m \^ r ( \ , m K* . 



ym B 
22/ 



(Cj.("iB>y) + Ci("iB,!/)) 



<i(m B ,y) ( 2 (m B ,y) 



ym B 



m B + m K * 

(1 + H^l)(£ ± ( mB;2/ ) _ _Kl^( mBi yfi 
m B y 

+ — )C±(m*,y) 



m B 



y 



C±(m B ,y) + C,i{m B ,y) 

m B 

q 2 



m B — m 



)Ci(™b, !/) + (! + 



2 ~) Ci(m B ,y) 

m B — m K * m B 



T 3 {q 2 ) = (±(m B ,y) - (1 - —{C\\{m B ,y) +(2(m B ,y)) - —(i(m B ,y) 



rn 



"IB 



m B 



5.1 



5.2 



5.3 



5.4 



5.5 



5.6 



5.7 

5.8 
5.9 

(5.10 

(5.11 
(5.12 
(5.13 
(5.14 

(5.15 

(5.16 
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Comparing the above forms with the ones given by eqs.(104)-(113) in Ref. jD]], it is ob- 
vious that the three functions Ca(q 2 ), (i(q 2 ) an d (2(q 2 ) should vanish, i.e., Ca(q 2 ) = Ci(<? 2 ) — 
C2(<? 2 ) — 0, in order to reproduce the results in LEET. For a quantitative comparison, we 
plot in Fig. 21 three curves (\\{q 2 ), (±(q 2 ) an d C(<? 2 ) as functions of momentum transfer q 2 , 
and in Fig.22 the three ratios (2(q 2 ) / (\\(q 2 ) , (i(q 2 ) / (±(q 2 ) and Ca(<? 2 )/C(<? 2 ) as functions of 
q 2 . It is seen that the three functions Ca(? 2 ), Ci(<? 2 ) an d (2(q 2 ) do have small but sizable 
contributions relative to the three nonzero functions (\\(q 2 ), (±{q 2 ) and ({q 2 ) in LEET. Nu- 
merically, they are about 20%, 10% and 10% respectively when q 2 = 0. For the ratios 
C2(<? 2 )/C||(<? 2 ) an d (A(q 2 )/((q 2 ), they are almost independent of the momentum transfer q 2 . 
For the ratio Ci(<? 2 )/Cl(<? 2 )) it ranges from 10% to 40% as q 2 increases. 

To be more clear, we directly plot in Figs. 23-29 the ratios between the form factors, i.e., 
FiLEET^q 2 ) I Fi(q 2 ) at center value of sq with T = 2.0GeV. Here F iEEE T{q 2 ) denote the form 
factors obtained from HQET/LEET relations, namely C4(<? 2 ) = Ci (<7 2 ) — ^(q 2 ) = 0, and 
Fi(q 2 )s are the form factors obtained in present paper, namely the three functions Ca(q 2 )) 
Ci(q 2 ) a nd (2(q 2 ) axe not zero and given in eq. (5.2), (5.5) and (5.6). It is seen that, except 
the ratio A leet/A w 0.78, the other ratios near large recoil point (q 2 — > 0) have the 
value F iLEET /Fi = 0.9 ~ 1.0. More explicitly speaking Uleet/I+ ~ Itleet/ fr ~ 0.9, 
A-lleet/Ai « 0.92, T lLEET « T 2LEET /T 2 « 0.96 and T 3LEET /T 3 « 0.98. So we can say 
the form factors F iLEET have the accuracy beter than 90% as comparation with F, except 
Aqleet whose accuracy is only about 78%. This exception can be explained as that the third 
term on the R.H.S. of the equation (5.10) contains the factor y/rriK*, which undoubtedly 
enhances several times the contribution of (i{q 2 ) when q 2 = 0. So when B decays to more 
light final mesons, say p, the LEET relation for Ao(q 2 ) will work worse. It is particularly 
noted that ratios J+leetI f+, Ti LEET /T 1 and T 3LEET /T 3 , which take value among the area 
0.85 ~ 0.90, 0.91 ~ 0.95 and 1.0 ~ 1.05 respectively in the whole q 2 , are almost independent 
of q 2 . The reason is clear since the ratios Ca(9 2 )/C(? 2 ) i n (5-7) and (2{q 2 ) / (\\(q 2 ) in (5.16) are 
almost independent of the momentum transfer q 2 and the coefficient of the wave function 
Ci( m B, y) is associated with the factors y/rriB in (5.14) and (5.16), which have a tendency to 
suppress the contribution of Ci{ m B, y) a s q 2 increases. While the other ratios, especially the 
ratio A 0LEET /A , get worse results for large q 2 , say q 2 ~ lbGeV 2 , which can be considered 
as far from large recoil. For instance Aq EEE t is only about half of A at q 2 ~ 15GeV 2 . 
This means that LEET becomes not appropriate for large q 2 region. The reason lies in 
that on one hand (i{q 2 ) / (±{q 2 ) counts notablly when q 2 is large. On the other hand, in 
the general form factor formulations (5.8), (5.10) and (5.15), there is a factor q 2 to enhance 
the contributions of Ca(q ,2 ) ! (2(q 2 ) an d Ci(<? 2 ) when q 2 increases. In addition, as can be seen 
from (5.9), (5.12) and (5.13), /tleet, A 2 leet an d Vleet are same as the corresponding 
one obtained by HQEFT in present paper. 

Besides the relations (3.13-3.16), LEET leads to three additional relations |L7| 

f T (q 2 ) = (l + ^)f + (q 2 ) (5.17) 
m B 

M f) = (1 - ^1.)—^ _ Mf)) + -K-i<ns + mK . ) 

ym B m K *(m B + m K *> Zym B 

To check this relations, we plot in Figs. 30-32 the three form factors fx, A and A\ obtained 
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from the direct calculations (i.e., left-hand side (LHS) ) and from the relations (i.e., right- 
hand side (RHS)). The numerical comparison of the LHS and RHS of the above equations 
shows that near the large recoil point (q 2 — > 0) the relations (5.17) and (5.19) both hold 
within 90% accuracy, but the relations (5.18) only hold about 80% accuracy. In particular, 
the deviations of LHS and RHS in relations (5.17) and (5.19) do not exceed 20% through 
whole q 2 , while that of relation (5.18) at best hold with 40%. 



VI. RESULTS FOR BRANCHING RATIOS 



The relevant branching ratios are able to be calculated with the form factors given above. 
The relevant decay width formulae have the following general forms fl] 



dT GWm% 



ds 



B \v;M 2 m{(\A'(s)\ 2 + \c'(s)\ 2 )(\(s) 



2 n 7T 5 



u(s) 



+ \C'{s)\ 2 4m 2 {2 + 2m 2 K - s) + Re{C'{s)D'{s)*)8m 2 {l - m 2 K ) + \D'{s)\Hm 2 s} (6.1) 
for B — ► Kll, and 

\Ms)\ 2 



dT G 2 F a 2 m% 2 

\V u Vib\ u(s 



m l\ , I 77V~\|2- M ' S ' 



ds 



2 n 7T 5 
1 

Am 2 



? -sX(s)(l + 2^) + \E(s)\ 2 s ? 



K* 



4m 2 K , 



\B(s)\ 2 (X(s) 
\C(s)\ 2 (m 



+8mUs + 2mf)) + \F(s)\ 2 (m ,, 



+ 8rh 2 K *(s - Arhf)) 



, + \G(s)\ 2 (X(s) -=^- + 4m/(2 + 2m 2 K , - g)) 



2m 2 K , 



Re{B{s)C{s)*){X{s) 



u[s) 



+Re(F(s)G(s)*)((\(s) 



U(S 



-){l-m 2 K ,-s)+Am 2 X{s)) 



-2^-X(s) \Re(F(s)H(sY) - Re(G(s)H(sY)(l - m 2 ) + ^-sX(s)\H(s) 
mr K * L J mr K * 



mt 



(6.2) 



for B — > K*ll, as well as 
L 



G 2 w am 2 rn? B 



\v:M 2 \c e r ff m(omi-m K , 



32tt 4 



(6.3) 



for B -> K* 7 . 

The functions appearing in the above decay width formulae are defined as [[TlflTI 



A'(s)=C e /f(s)f + (s) + 

B'(s)=C e / f (s)f4s) 
C'(s)=C 10 f + (s), 



-C?'f T {8) 



1 + m K 



2rhb 



(1 - m K )C e 7 ff f T (s) 



(6.4) 

(6.5) 
(6.6) 
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D'(s)=C 10 f-(s), (6.7) 

A(s) = —^—C e 9 ff (s)V(s) + ^C^T^s), (6.8) 

B(s) = (1 +m K .)[C% ff (s)Ai(s) + ^(l-m K *)C e 7 ff T 2 (s)}, (6.9) 

C(s) = \ 2 [(1 - m K *)Cl ff (s)A 2 (s) + 2m b C e 7 ff (T 3 (s) + 1 ~ ?*• T 2 (s))], (6.10) 
1 — * s 

D(s) = l[C e 9 ff {s)({l + m K .)A x (s) - (1 - m*.)^) - 2m K *A (s)) - 2m b C e / %(§)}, (6.11) 
s 

E(s) = I Cio^(s), (6.12) 
1 + mjc* 

= (l+m K .)C 1 oA 1 (§), (6.13) 

G(s) = \ C 10 A 2 (s), (6.14) 
1 + m K * 

H{§) = |Ci (s)[(l + rh K ,)A 1 (s) - (1 - m K *)A 2 {s) - 2m K *A {s)}, (6.15) 

A(s) = l+m^^ + s 2 -2s-2m^^(l + s), (6.16) 



6(S)=yA(S)(l-4^), (6.17) 
C 9 e// (s) = C 9 + ff (m c , s)(3d + C 2 + 3C7 3 + C 4 + 3C7 5 + C 6 ) - ^(1, S)(4C 3 + 4C7 4 + 3C 5 + C 6 ) 

- ^(0, s)(C 3 + 3C 4 ) + ^(3C7 3 + C 4 + 3C7 5 + C 6 ), (6.18) 

8 8 m b 4 , 4 . 
S (0, S ) = -- 5 l n -- 5 l 11S + -», (6.19) 

, „, 8 8, m b 8, 4 

g (z, s = m in z H — x 

yv ' ; 27 9 /i 9 9 

9 v 2 arctan , = , for x = =+ > 1. 

with s = -^.vtlkk* = mK ' K * rhb c = and rhi = — . For the Wilson coefficients Cj, we 
take the results calculated in the naive dimensional regularization (NRD) scheme [42] and 
their values are listed in Table 2. 



Table 2. Values of the Wilson coefficients with choosing the renormal- 

y 7 



ization scale at fj, = m b = 4.8GeV. Here, G% * = Cj — C5/3 — C%. 



Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


CV 


C 7 e// 


c 9 


C10 


-0.248 


1.107 


0.011 


-0.026 


0.007 


-0.031 


-0.342 


-0.313 


4.344 


-4.669 



It is seen that once the decay rates and the form factors are precisely determined, one is 
able to extract the CKM matrix elements. Here we may use the current reasonable values 
of CKM matrix elements extracted from other processes and unitarity of the CKM matrix 
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to predict the branching ratios for the B meson rare decays. When taking the CKM matrix 
element IV^V^I = 0.0385, we present our results in Table 3. For comparison, we also list the 
results given in |lj and the ones from experimental measurements. It is seen that within the 
uncertainties, our prediction is consistent with the experimental results. 

Table 3. Branching ratios (Br) for B rare decays in standard model 
(SM). In deriving the branching ratios we have used the lifetime of B me- 
son: tb = 1.65ps. Where the errors mainly come from the uncertainties 
of the threshold energy sq, the Borel parameter T and the parameters 
in (4.3). 





present values 


values in |j] 


experiment 


B -> Kee 
B — > Kfijl 


(0.84i°;^)xl0- 6 


(0.57i°;^)xl0- 6 


(0.75l[J;|±0.09)xlO- 6 |43| 


B -> Ktt 


(1.73--) x 10-' 


(1.3^;?) xlO" 7 




B -> K*-y 


(5.47i|^)xl0- 5 


(5.37iS)xlO" 5 A 


(4.551&6B ± 0.34) xlO" 5 AA |44 
(4.96± 0.67±0.45)xl0~ 5 AA [45 
(5.7 ± 3.3)xl0" 5 H 


B -> K*ee 


(1.86i^)xl0- 6 


(2.3i^)xl0- 6 


(2.08^i;^37) x io-° m 


B -> K*fifl 


(1.78+^) x 10~ e 


(1.9^) xlO- 




B K*rf 


(1.68i^)xl0- 7 


(1.9i^)xl0- 7 





A The value is obtained through (6.3) by using the form factor Ti(0) = 
0.379+g;g|| presented in Ref. Q. 
AA : The value is for B° -» K°* 7 . 



VII. SUMMARY 

In summary, we have studied the B meson rare decays in the framework of HQEFT. The 
semileptonic type and penguin type form factors have been derived by using LCSR method 
in HQEFT. It has been seen that the heavy quark expansion brings a much simplification to 
B meson rare decays. Isgur-Wise relations among the semileptonic type and penguin type 
form factors have been proved to hold for the whole momentum transfer region at the leading 
order I/tuq expansion. As a consequence, all the form factors can be neatly characterized by 
a set of wave functions (A, B and Lj(z = 1, 2, 3, 4)) at the leading order of 1/tuq expansion. 
Furthermore from our quantitative discussion, it is obvious that LEET is a valid method for 
heavy-to-light transition. LEET relations hold within 80% accuracy at large recoil point on 
the whole, and most of them even hold better than 90% accuracy. Moreover, our numerical 
prediction is consistent with the experimental results. We then conclude that the branching 
ratios of B meson rare decays can be reasonably predicted based on LCSR approach within 
the framework of HQEFT. Nevertheless, in order to match the expected measurements of B 
factories in the near future, a more accurate calculation of the form factors for the B meson 
rare decays is urgently required. 
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FIGURES 



, , , , , , T [GeV] 

1 1.5 2 2.5 3 3.5 

Fig.l The form factors /+ as functions of 
the Borel parameter T for different values 
of the threshold sq. The dotted, solid and 
dashed curves correspond to sq = 2.1,2.4 
and 2.7GeV respectively. The momentum 
transfer is q 2 = 0. 
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Fig. 2 Same as Fig.l but for /_. 
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Fig.3 Same as Fig.l but for fx- 
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Fig. 4 The form factors Aq as functions of 
the Borel parameter T for different values 
of the threshold so- The dotted, solid and 
dashed curves correspond to so = 1.8,2.1 
and 2.4GeV respectively. The momentum 
transfer is q 2 = 0. 
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Fig. 5 Same as Fig. 4 but for A\. 



Fig. 6 Same as Fig. 4 but for A2. 
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Fig.7 Same as Fig. 4 but for V. 



Fig. 8 Same as Fig. 4 but for T\. 
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Fig. 9 Same as Fig. 4 but for T<i- 
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Fig. 10 Same as Fig. 4 but for T3. 
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Fig. 11 The form factors /+ as functions 
of q 2 for different values of the threshold 

sq. The dotted, solid and dashed curves Fig.12 Same as Fig.ll but for /_. 

correspond to so = 2.1,2.4 and 2.7GeV 

respectively. 
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Fig.13 Same as Fig.ll but for fx- 



Ai 




Fig. 15 Same as Fig. 14 but for A\. 




Fig. 17 Same as Fig. 14 but for V. 
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Fig. 14 The form factors Aq as functions 
of q 2 for different values of the threshold 
so- The dotted, solid and dashed curves 
correspond to sq = 1.8,2.1 and 2.4GeV 
respectively. 
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Fig. 16 Same as Fig. 14 but for A2- 




Fig. 18 Same as Fig. 14 but for Ti. 
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Fig. 19 Same as Fig. 14 but for T 2 . 



Fig.20 Same as Fig. 14 but for T 3 . 
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Fig.21 C||(? 2 ), C±(q 2 ) and ((q 2 ) at 
center value of so with T = 2.0GeV. 
The dotted, dashed and solid curves 
correspond to C||(</ 2 )> C±(l 2 ) and 
C(q 2 ) respectively. 
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Fig.22 C2(9 2 )/C||(S 2 ), Ci(q 2 )/C±(q 2 ) and 
(,a{<1 2 ) / (,{q 2 ) at center value of sq with 
T = 2.0GeV. The dotted, dashed and 
solid curves correspond to (,2{q 2 ) / C,\\{q 2 ) i 
(i{q 2 )/(±(q 2 ) and ( A (q 2 ) / 'CO? 2 ) respectively. 
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Fig.23 f+LEET(q 2 )/f+(q 2 ) at 
s = 2.4GeV with T = 
2.0GeV. 
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Fig. 24 Same as Fig.23 but for 
f~LEET(q 2 )/f-{q 2 ). 
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Fig.25 A l££t(9 2 )Mo(9 2 ) at 
s = 2.1GeV with T = 
2.0GeV. 



Fig. 26 Same as Fig.25 but for 
A\leet (Q 2 )/A\(q 2 ). 
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Fig. 27 Same as Fig.25 but for 

TiLEETiq 2 )/^ 2 ). 



Fig. 28 Same as Fig.25 but for 

T 2L EET{q 2 )/T 2 (q 2 ). 
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Fig. 29 Same as Fig.25 but for 
T 3LE ET{q 2 )/T 3 {q 2 ). 
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Fig.30 LEET relation (5.17) for 
B -» K at s = 2.4GeV and 
T = 2.0GeV. The dashed and solid 
curves correspond to LHS and RHS 
of the equation respectively. 
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Fig.31 LEET relation (5.18) for 
B — > K* at so = 2.1GeV and 
T = 2.0GeV. The dashed and solid 
curves correspond to LHS and RHS 
of the equation respectively. 
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Fig.32 LEET relation (5.19) for 
B K* at s = 2.1GeV and 
T = 2.0GeV. The dashed and solid 
curves correspond to LHS and RHS 
of the equation respectively. 
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